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CHAPTER 1. INTRODUCTION 


1.1.3 Sobolev spaces 

In this section we introduce a class of spaces, called Sobolev spaces (after the Russian 
mathematician S.L. Sobolev), which play an important role in modern differential 
equation theory. Before we give the precise definition of a Sobolev space, we intro¬ 
duce the concept of weak derivative. 

Suppose that u is a smooth function, say u G C fc (12), with 12 an open subset of 
M n , and let v G Cp 0 ^); then the following integration-by-parts formula holds: 

/ D a u(x) ■ v(x) dx = (—1)'"' f u(x) ■ D a v(x) dx, \a\ < k, 

Jn Jn 

\/v G (75° (fi). 

Note that all terms involving integrals over the boundary of 12, which arise in the 
course of integrating by parts, have disappeared because v and all of its derivatives 
are identically zero on the boundary of 12. This identity represents the starting point 
for defining the concept of weak derivative. 

Now suppose that u is a locally integrable function defined on 12 (i.e. u G Lficu) 
for each bounded open set oj, with a) C 12). Suppose also that there exists a function 
w a , locally integrable on 12 and such that 


f w a (x) ■ v(x) dx = (— l)l“l f u(x) ■ D a v(x) dx Vr G C 0 °°(l 2); 

Jn Jn 

then we say that w a is a weak derivative of the function u of order |a| = aq +... + 
a n , and we write w a = D a u. In order to see that this definition is correct it has to 
be shown that if a locally integrable function has a weak derivative then this must be 
unique; we remark that this is a straightforward consequence of DuBois Reymond’s 
lemma 5 . Clearly, if u is a sufficiently smooth function, say u G C k {1 2), then its weak 
derivative D a u of order \a\ < k coincides with the corresponding partial derivative 
in the classical pointwise sense, 

d l“l-u 

(9a;" 1 ... <9.x'" ri ’ 

In order to simplify the notation, we shall use the letter D to denote classical as 
well as weak derivatives; it will always be clear from the context (by considering the 
smoothness of the function differentiated) which of the two is implied. 


Example 5 Let 12 = M 1 , and suppose that we wish to determine the weak first 
derivative of the function u(x) = (1 — |x|) + defined on 12. Clearly u is not differen¬ 
tiable at the points 0 and ±1. However, because u is locally integrable on 12, it may, 


°DuBois Reymond’s lemma: Suppose that w is a locally integrable function defined on an 
open set Cl, Q C R n . If 


[ w(x)v(x) dx = 0 for all v in Cq°(CI) 
Jn 


then w(x) = 0 for almost every x € Cl. 
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nevertheless, have a weak derivative. Indeed, for any v G Cg°(Q), 


r+oo 


r+oo 


u{x)v'(x) dx — (1 — |a;|)_)_t/(a;) dx — (1 — |a;|)t/(a;) dx 


'-i 


= / (1 + x)v\x) dx + / (1 — x)v'(x) dx 


'-i 


= — I v(x) dx + (1 + x)v(a;)| 0 x + j v(x) dx + (1 — x)v(x)\]. =0 

w(x)v(x) dx , 


J -1 JO 

r*0 pi /»+oo 


= J (—l)v(x)dx + J \-v(x)dx = 


where 


w(x) = 


0, x < —1, 

1, are (-1,0), 

-1, are (0,1), 

0, x > 1. 


TTras, the piecewise constant function w is the first (weak) derivative of the contin¬ 
uous piecewise linear function u, i.e. w — vl = Du. o 

Now we are ready to give a precise definition of a Sobolev space. Let k be a 
non-negative integer and suppose that p G [1, oo]. We define (with D° denoting a 
weak derivative of order |a| ) 

W£(Q) = {«G Lp(ft) : G L p (0), |a| < A:}. 

Wp(Q) is called a Sobolev space of order k ; it is equipped with the (Sobolev) norm 


i/p 


PlIwKn) 


|a|<fc 






when 1 < p < oo 


and 


Mlw* («) : = J] H^^IUocin) when p = oo. 

|a:|<fc 


Letting, 


i/p 




“„iip 

L v (£l) 


\a\=k 


for p G [1, oo), we can write 


i/p 




Ww£(Q) ~ ( X 




,1=0 
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Similarly, letting 


we have that 


u \ W*,(Sl) 


Y. ii jC,aM iUoo(n), 

\a\=k 


k 

u \\w^(n) = Y, \ u \wLO)- 

3=0 


When k > 1, |*|w^(si) is called the Sobolev semi-norm 6 on VF* (ft). 

An important special case corresponds to taking p — 2; the space W^(ft) is then 
a Hilbert space with the inner product 


(u,v) W k {n) := Y, (D a u, D a v). 

|a|<fc 


For this reason, we shall usually write H k (D,) instead of W 2 (Pl). 

Throughout these notes we shall frequently refer to the Hilbertian Sobolev spaces 
if 1 (ft) and H 2 {Vt). Our definitions of (ft) and its norm and seminorm, for p = 2, 
k — 1, give: 


if 1 (ft) = 

Mrrqn) = 

Similarly, for p = 2 and k 
H 2 {0) = 


u e L 2 (ft) : w— e L 2 (ft), J = 


dx 


u 


+ E 


du 


L2O) -r H q x ,Wl 2 {Q) 


3 =1 


} 


1/2 


Y 11 du " 2 


.3 =1 

= 2, 


1 


1/2 


9x t - ll£a ^ 


du 


u e L 2 (0) : — e L 2 (0), j = 1,..,, n, 


d 2 u 

dxidxj 


G L 2 (ft), i,j = , 


«||rr 2 (n) - {IkllL(a) + XI H^rllL(n) 

j=i 3 

d 2 u 
dxidxj 


6 When k > 1, | • \w£(n) is on ly a semi-norm rather than a norm because if |u| w£(n) = 0 for 
u £ Wp(Cl) it does not necessarily follow that u(x) = 0 for almost every x in ft (all that is known 
is that D a u(x ) = 0 for almost every x € f2, |a| = k), so | • \w k (Q) does not satisfy the first axiom 
of norm. 


E 

t,j =1 
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u\h 2 (Q.) 



d 2 u 2 
dxidxj 


1/2 


Finally, we define the special Sobolev space Hq(Q) as the closure of Cq°(Q) in the 
norm of || • Hnpo); in other words, Hq(Q) is the set of all u G iF 1 (h2) such that u 
is the limit in H 1 ^Ll) of a sequence {u m }™ =1 with u m G Cq° (hi). It can be shown 
(assuming that dLl is sufficiently smooth) that 

H ( ] (O) = {mG iF 1 (Q) : u — 0 on <9f2}; 

i.e. Hq(Q) is, in fact, the set of all functions u in i/ 1 (h2) such that u = 0 on dLl, the 
boundary of the set hi. We shall use this space when considering a partial differential 
equation that is coupled with a homogeneous (Dirichlet) boundary condition: u — 0 
on dLl. We note here that Hq(Q.) is also a Hilbert space, with the same norm and 
inner product as ff 1 (f2). 

We conclude the section with the following useful result. 


Lemma 2 (Poincare-Friedrichs inequality) Suppose that Q is a bounded open set in 
M n (with a sufficiently smooth boundary 1 dQ) and let u G Hq (Q); then there exists 
a constant c*(fi), independent of u, such that 



u(x) | 2 dx < c* 



( 1 . 2 ) 


Proof As any function u G H ( \ (Q) is the limit in if 1 (H) of a sequence {tt m }“ =1 C Cq°(Q), 
it is sufficient to prove this inequality for u G C'o°(H). 

In fact, to simplify matters, we shall restrict ourselves to considering the special case 
of a rectangular domain = (a, b) x (c, d) in M 2 . The proof for general is analogous. 
Evidently 


u(x, y) 


r x 

u ( a ,y) +J 



du 

dx 


(f,y) dC, 


c < y < d. 


Thence, by the Cauchy-Schwarz inequality, 


| u(x, y)\ 2 dx d y 


= fflf^ymdx 

< J\x-a)dx^ j\^{Cy)\ 2 dfdy^j 
= \{b~ a) 2 j^\^{x,y)\ 2 dxdy. 


7 Say, Pis a polygonal domain in R 2 or a polyhedron in R 3 . 
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Analogously, 


[ \u(x,y)\ 2 dxdy < \{d - c) 2 [ (x, y)\ 2 dx dy. 

Jn 2 J n dy 


By adding the two inequalities, we obtain 


\u(x,y)\ 2 dxdy < c*^ dxd V^ 


where c* = ((^ + . ■ 

For further reference, we note that if O = (0, l ) 2 C M 2 then c* = similarly, if 
O = (0,l) CR then c* = 


1.2 Weak solutions to elliptic problems 

In the first part of this lecture course we shall focus on boundary value problems for 
elliptic partial differential equations. Elliptic equations are typified by the Laplace 
equation 


A u = 0, 

and its non-homogeneous counterpart, Poisson’s equation 


-Aw = /, 


where we used the notation 

71 

1=1 L 

for the Laplace operator. 

More generally, let 17 be a bounded open set in M n , and consider the linear 
second-order partial differential equation 

^^dx-V^^dx) + ^Ux)— + c{x)u = f{x), xen > ( L3 ) 

where the coefficients a^-, &*, c and / satisfy the following conditions: 

a»j e C 1 ^), i,j = l,...,n- 
bi e C(f2), i — 1 ,..., n; 
ceC(fi), /eC(fi), 

n n 

£ n )e®", I6fi; 


and 


(1.4) 
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here c is a positive constant independent of x and £. The condition (1.4) is usually 
referred to as uniform ellipticity and (1.3) is called an elliptic equation. 

In problems that arise in applications equation (1.3) is usually supplemented by 
one of the following boundary conditions, with g denoting a given function defined 
on <90: 


(a) 

(b) 


(d) 


u = g on <90 (Dirichlet boundary condition); 


= g on <90, where v denotes the unit outward normal vector to <90 (Neu¬ 
mann boundary condition); 

+ au — g on <90, where a(x ) > 0 on (90 (Robin boundary condition); 

A generalisation of the boundary conditions (b) and (c) is 

<9 u 


£ 

1,3= 1 


dij— — cos ctj + cr(x)u = g on <90, 

t/X j 


where oij is the angle between the unit outward normal vector v to <90 and 
the Xj axis (Oblique derivative boundary condition). 


In many physical problems more than one type of boundary condition is imposed 
on dLl (e.g. <90 is the union of two disjoint subsets <90 x and <90 2 , with a Dirichlet 
boundary condition on <90 1 and Neumann boundary condition on <902). The study 
of such mixed boundary value problems will not be pursued in these notes. 

We begin by considering the homogeneous Dirichlet boundary value problem 


v - ^ <9 f du \ . . du . . , . 

2^ ( r/ 'R Ar . ) +2^b i {x)— + c{x)u = f{x), xeo, 

i,j= 1 3 \ 1 / j=l * 


u = 0 on <90, 


(1.5) 

( 1 . 6 ) 


where a l3 , bi, c and / are as in (1.4). 

A function u G C 2 (0) flC'(O) satisfying (1.5) and (1.6) is called a classical solu¬ 
tion of this problem. The theory of partial differential equations tells us that (1.5), 
(1.6) has a unique classical solution, provided that a^, bi, c, / and <90 are sufficiently 
smooth. However, in many applications one has to consider equations where these 
smoothness requirements are violated, and for such problems the classical theory is 
inappropriate. Take, for example, Poisson’s equation with zero Dirichlet boundary 
condition on O = (—1, l) n in 

-Au = sgn (| - \x\) , x G O, 1 . , 

u — 0, x G <90. J 

This problem does not have a classical solution, u G C 2 (0) 51(0(0), for otherwise Au 
would be a continuous function on O, which is not possible because sgn(l/2 — |x|) 
is not continuous on O. 
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In order to overcome the limitations of the classical theory and to be able to 
deal with partial differential equations with “non-smooth” data, we generalise the 
notion of solution by weakening the differentiability requirements on u. 

To begin, let us suppose that u is a classical solution of (1.5), (1.6). Then, for 
any v G C'o(fi), 




■ v dx 


c(x)uv da; 


+ 


n p 

e h 

i =i 


du 


x 


dxi 
f(x)v(x) da;. 




■ v dx 


Upon integration by parts in the first integral and noting that v — 0 on <90, we 
obtain: 




du dv 
dx.i dxj 


da; 


c(x)uv da; 


E 

i= 1 




f(x)v(x) da; 


v dx 

Vu G (7o(n). 


In order for this equality to make sense we no longer need to assume that u G U7 2 (0) : 
it is sufficient that u G L 2 (Q) and du/dxi G L 2 (0), i — 1,..., n. Thus, remembering 
that u has to satisfy a zero Dirichlet boundary condition, it is natural to seek u in 
the space H ( \ (Q), where, as in Section 1.1.3, 


Hq(Q) = {mG L 2 (D) : —— G L 2 (fl), i = 1, ..., n, u = 0 on <9h2}. 

C J ih n 


Therefore, we consider the following problem: find u in Hq(Q) such that 



du dv 


®ij (*^) o 


dxi dx , 


n p 

da; + V / bi 

i =i 


(x) w—v dx 

U X j 


+ / c(x)uvdx— / f(x)v(x)dx Vr G Cq( Q). 


(1.7) 


We note that Cq (0) C H ( \ (fl), and it is easily seen that when u G H ( \ (Q) and 
v G Hq (O), (instead of v G C { ) (O)), the expressions on the left- and right-hand side 
of (1.7) are still meaningful (in fact, we shall prove this below) 8 . This motivates the 
following definition. 


8 Note further that since the coefficients a, : j no longer appear under derivative signs in (1.7), 
it is not necessary to assume that € C ,1 (fi); a,y G LoUfl) will be seen to be sufficient. Also, 
the smoothness requirements imposed on the coefficients 6,; and c can be relaxed: bt G Aoo (ff) for 
i = 1,..., n and c G Aoo(f2) will suffice. 
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Definition 1 Let aij G L 00 (Q) f i,j = 1 ,...,n, bi G L 00 (ll) 7 i = 1 c G 

Loo (fl), and let f G L- 2 (Ll). A function u G Hq(Q) satisfying 




du dv 


dx 


dxi dxj 
c(x)uv dx 




f(x)v(x) dx 


v dx 

Vu G Hq(Q) 


( 1 . 8 ) 


Jn Jn 

is called a weak solution of (1-5), (1.6). All partial derivatives in (1-8) should be 
understood as weak derivatives. 


Clearly if u is a classical solution of (1.5), (1.6), then it is also a weak solution 
of (1.5), (1.6). However, the converse is not true. If (1.5), (1.6) has a weak solution, 
this may not be smooth enough to be a classical solution. Indeed, we shall prove 
below that the boundary value problem (*) has a unique weak solution u G Hq(Q), 
despite the fact that it has no classical solution. Before considering this particular 
boundary value problem, we look at the wider issue of existence of a unique weak 
solution to the more general problem (1.5), (1.6). 

For the sake of simplicity, we adopt the following notation: 



and 


l(v) 


f(x)v(x) dx. 


With this new notation, problem (1.8) can be written as follows: 


find u G such that a(u,v ) = l(v) Vv G i/g(fl). 


( 1 . 10 ) 


( 1 . 11 ) 


We shall prove the existence of a unique solution to this problem by exploiting the 
following abstract result from Functional Analysis. 

Theorem 1 (Lax & Milgram theorem) Suppose that V is a real Hilbert space equipped 
with norm || • ||y. Let a(-, •) be a bilinear functional on V x V such that: 

(a) 3c 0 >0 Vv G V a(v,v) > c 0 ||u||y, 

(b) 3ci > 0 \/v,W G V |a(iu, l?)| < Ci||w||y||u||y, 
and let /(•) be a linear functional on V such that 
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Then, there exists a unique u E V such that 

a(u,v)=l(v) WvEV. 

For a proof of this result the interested reader is referred to the books: P. Ciarlet: 
The Finite Element Method for Elliptic Problems, North-Holland, 1978; K. Yosida: 
Functional Analysis, Reprint of the 6th ed., Springer-Verlag, 1995. 

We apply the Lax-Milgram theorem with V = Hq(Q) and || • ||y = || • ||#i(Q) 
to show the existence of a unique weak solution to (1.5), (1.6) (or, equivalently, to 
(1.11)). Let us recall from Section 1.1.3 that Hq(Q.) is a Hilbert space with the inner 
product 


{w,t>) H i(n)= / wvdx + 


in 


n 

E 

1=1 


dw dv 
dx/dx, X 


-j /n 

and the associated norm ||iu||ffi(n) = (w, w)^ 1(T2 ). Next we show that a(-, •) and /(•), 
dehned by (1.9) and (1.10), satisfy the hypotheses (a), (b), (c) of the Lax-Milgram 
theorem. 

We begin with (c). The mapping v l(v) is linear: indeed, for any a,/3 G M, 


l(av i+/3v 2 ) = / f(x)(av 1 (x) + (3v 2 {x)) dx 


— a f(x)vi(x)dx + /3 / f(x)v 2 (x)dx 

Jn Jn 

= al(vi) + 0l(v 2 ), Vi, v 2 G Hq(D); 


so /(•) is a linear functional on Hq (Q). Also, by the Cauchy-Schwarz inequality, 
\l(v)\ = | f f(x)v(x) dx\ < f f \f(x)\ 2 dx^j f f \v(x)\ 2 dx^j 

Jn \Jn J \Jn J 

= ll/IU 2 (n)|M|L 2 (Q) < ||/||L 2 (n)lkl|ifi(n), 


for all v G Hq(Q), where we have used the obvious inequality |M|z, 2 (n) < ||^||jy 1 (V 2 )- 
Letting c 2 = ||/||z, 2 (n), we obtain the required bound. 

Next we verify (b). For any fixed w G H ( \ (Q), the mapping v i—> a{y, w ) is linear. 
Similarly, for any hxed v G Hq(Q), the mapping w i-G a{v,w) is linear. Hence a(-, ■) 
is a bilinear functional on H ( \ (H) x Hq(Q). Applying the Cauchy-Schwarz inequality, 
we deduce that 


\a(w,v)\ < y y max |a^-(x)| 


i,j= 1 


dw dv 
2 dxi dxj 

Tl /> O 

+ ^ max | bj (x) | | j r dx\ 


dx | 


i =1 




in 
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+ max|c(x)|| / w(x)v(x)dx\ 

J n 


* He 

\*j=i 

n 

+E 




i —1 


n dx 

I \^\ 2 dxY 2 

' n QXi ) 


L Kdx )' 


f \ ( f \ ^ b2 

I \w \ 2 dx\ (J \v \ 2 dx\ 


< d 


{(bM’-±XL 

{(IM’AM. 


dw \2, y /2 

n dxi X ) 


^,2 d y /2 
^ / 


where 


c = max < max max \aij(x)\, max max |6j(x)|, max |c(x)| > . 

l<i,j<n xefl 1 <i<n X £fl x£fl 


By further majorisation of the right-hand side in (1.12) we deduce that 

1/2 

I /' . - /• rhn _ I 

| a(w, v)| < 2 nd ■ 


x 


{ [ \w\ 2 dx + Y] [ 
[Jo i=1 Jn 


9W \2 A \ 
a ^ J 
1/2 


dv 

dxj 


dx 


so that, by letting C\ = 2 nd, we obtain inequality (b): 

I a(w,v)\ < ci||w||^i (r2) ||u||^i (n) . 


( 1 . 12 ) 


(1.13) 


It remains to establish (a). To do so, we shall slightly strengthen the smoothness 
requirements on the coefficients b t by demanding that bi G W^II) (see, however, 
Remark 4 at the end of this chapter). Using (1.4), we deduce that 


Tl n r\ TL n 

a{v,v)>c^ / \g^\ 2 dx + J2 / b > 

i=i J n ° Xi i= i 


(x)}--^-(v 2 ) dx + f c(x)\v\ 2 dx, 

2 oXi Jq 


where we wrote • v as §//“~H ;2 )- Integrating by parts in the second term on the 
right, we obtain 
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Suppose that i — 1,..., n, and c satisfy the inequality 


Then 


c(x) 



2—1 


dbj 

dxi 


> 0 , 


x G fb 



(1.14) 


(1.15) 


By virtue of the Poincare-Friedrichs inequality stated in Lemma 1.2, the right-hand 
side can be further bounded below to obtain 


a(v,v ) > 


— / \v\ 2 dx. 
c * Jn 


Summing (1.15) and (1.16), 


a (v, v) > c 0 



(1.16) 


(1.17) 


where Co = c/(l + c*), and hence (a). Having checked all hypotheses of the Lax- 
Milgram theorem, we deduce the existence of a unique u G Hq(D) satisfying (1.11); 
consequently, problem (1.5), (1.6) has a unique weak solution. We encapsulate this 
result in the following theorem. 


Theorem 2 Suppose that aij G L^Ul), i,j = 1, ... ,n, b, G W^,(fi), i = 1, ... ,n, 
c G L 0O (n), / G L 2 (D), and assume that (1-4) and (1-14) hold; then the boundary 
value problem (1.5), (1-6) possesses a unique weak solution u G H) (Q). In addition, 

llwllnqo) < —||/||l 2 (o)- (l- 18 ) 

Co 

Proof We only have to show (1.18) as the rest of the theorem has been proved above. 
By (1.17), (1.11), the Cauchy-Schwarz inequality and recalling the definition of || • H^qn), 

c olMI/p(fi) < a(u,u) = l(u ) = (f,u) 

< |(/, «)| < ||/IU 2 (n)l(«llL 2 (Q) 

< ll/IU 2 (Q)ll w ll// 1 (n)- 

Hence the desired inequality. ■ 

Now we return to our earlier example (*) which has been shown to have no 
classical solution. However, applying the above theorem with atj(x) = 1, i = j, 
aij(x ) = 0, i ^ j, 1 < i,j < n, b,(x) = 0, c(x) = 0, f(x) = sgn(| — |x|), and 
Q = (—l,l) n , we see that (1.4) holds with c = 1 and (1.14) is trivially fulfilled. 
Thus (*) has a unique weak solution u G H) (H) by Theorem 2. Similar results 
are valid in the case of Neumann, Robin, and oblique derivative boundary value 
problems, as well as mixed problems. 
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Remark 2 Consider, for example, the following Dirichlet-Neumann mixed bound¬ 
ary value problem: 

f in 

0 ora Ti, 

9 on r 2 , 

where is a non-empty, relatively open subset of dO and r\ U T 2 = <9f2. We shall 
suppose that f G L 2 (L1) and that g G L 2 (T 2 ). Following a similar reasoning as in the 
case of the Dirichlet boundary value problem, we consider the special Sobolev space 


—A u = 
u = 
du 
dv 


= {v G if 1 (ft) : v = 0 ora TJ, 

and define the weak formulation of the mixed problem as follows: find u G Hq T (Q) 
such that 

a(u,v ) = l(v) for all v in Hq T (Q), 


where we put 


a(u, v) 



du dv 
dxidxl x 


and 


l(v) = / f(x)v(x)dx + / g(s)v(s) ds. 
Jn J r 2 


Applying the Lax-Milgram theorem with V = Hq V (LI), the existence and uniqueness 
of a weak solution to this mixed problem easily follows. o 


Remark 3 Theorem 2 implies that the weak formulation of the elliptic boundary 
value problem (1.5), (1-6) is well-posed in the sense of Hadamard; namely, for each 
f G L 2 (Q) there exists a unique (weak) solution u G (Q), and “small” changes in 
f give rise to “small” changes in the corresponding solution u. The latter property 
follows by noting that if U\ and u 2 are weak solutions in Hq(Q) of (1.5), (1.6) 
corresponding to right-hand sides fi and f 2 in L 2 (Q), respectively, then U\ — u 2 
is the weak solution in H((Q) of (1.5), (1.6) corresponding to the right-hand side 
/ 1 -/ 2 G L 2 (0). Thus, by virtue of (1.18), 

IK _ u 2 \\ H i(Q) < —\\fi - /2||z, 2 (n), (1-19) 

Co 

and hence the required continuous dependence of the solution of the boundary value 
problem on the right-hand side, o 

Remark 4 The requirement bi G Wf.)!!) in Theorem 2 can be relaxed to the original 
assumption bi G L 00 (H), i = 1,... ,n. To see this, note that the smoothness require¬ 
ments on bi are unrelated to the verification of condition (c) in the Lax-Milgram 



